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Abstract 

We investigate how strong a hypothetical bound state of two neutrons would affect different observables 
in the neutron-deuteron reactions. To that aim we extend our momentum space scheme of solving three- 
nucleon Faddeev equations to incorporate in addition to the deuteron also the dineutron bound state. 
We discuss effects induced by dineutron on the angular distribution of the neutron-deuteron elastic scattering 
and cross sections of the deuteron breakup. A comparison to the available data for neutron-deuteron total 
cross sections and elastic scattering angular distributions cannot decisively exclude a possibility that the 
two neutrons can form ^5*0 bound state. However, the strong modifications of a final-state- interaction peak 
of the neutron-deuteron breakup when changing from negative to positive values of the neutron-neutron 
scattering length seems to exclude existence of dineutron. 

PACS numbers: 21.45.-v, 21.45.Bc, 25.10.+S, 25.40. Dn 
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I. INTRODUCTION 



Investigation of the neutron-deuteron (nd) elastic scattering and the deuteron breakup reac- 
tion jl| revealed a number of discrepancies between data and theoretical predictions based on 
modern nucleon-nucleon potentials such as AV18 CD Bonn [3], and Nijml, 2 and 93 [J] or 
on the nuclear forces derived in the framework of chiral perturbation theory jsl- These potentials 
describe very accurately all existing nucleon-nucleon (NN) data as expressed by the value of 
per data point close to ~ 1. Some of these discrepancies can be explained when in adddition 
to pairwise interactions also three-nucleon forces (3NF's) are included in the 3N Hamiltonian. 
However, some persistently avoid explanation and they reveal high insensitivity to the underlying 
dynamics, especially to different choices among available 3NF's. The neutron-neutron (nn) quasi- 
free-scattering (QFS) configuration in the complete nd breakup is one such an example. Another 
is the symmetrical space-star (SST) geometry in that reaction. The nn QFS refers to the kine- 
matical configuration in which the outgoing proton is at rest in the laboratory system. In case 
of SST three outgoing nucleons have equal magnitudes of momenta, which in the three-nucleon 
cm. system form a plane perpendicular to the incoming nucleon momentum with the angle of 
120° between two consequitive momenta. In QFS and SST configurations theoretical predicitons 
drastically underestimate by about ~ 20% the data. This together with the fact that cross sec- 
tions in these configurations are dominated by the ^Sq and "^Si contributions lead to suspect that 
somethiiig is wrong with the ^Sq nn force and that maybe two neutrons can even form a bound 
state ^, Ql . 

That motivated us to investigate what consequences the existence of the dineutron would 
have on different observables in nd reactions and to what degree the available data for nd reactions 
allow for such a bound state. Also we would like to see if existence of dineutron could help in 
resolving the discrepancies present in QFS and SST nd breakup configurations. 

In section |TI] we extend our formulation of the momentum space treatment of the 3N Faddeev 
equations to include, in addition to the deuteron, also the ^5*0 bound state of two neutrons. Chang- 
ing the strength of ^Sq nn interaction of the CD Bonn potential we produce a number of forces 
which allow for two neutrons being bound with different dineutron binding energy. In section IIIII 
we present theoretical predictions based on solution of 3N Faddeev equations and compare them 
to the available nd data. We summarize and conclude in section ITVl 
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II. FADEEV EQUATIONS WITH DINEUTRON 



We shortly present the basics of our momentum space treatment of 3N Faddeev equations 



and calculations of the transition operators for different reactions in t 



re 3N continuum based on 
We put emphasis on 



solutions of these equations. For detailed presentation we refer to 
changes in the standard approach in the case when, beside the deuteron, also one additional bound 
state appears in some partial wave. 

For calculation of processes initiated from a state I'&i,! >= {qoi 4'd >, which describes the 
neutron moving with the relative momentum % with respect to the deuteron of the wave function 

one needs the state |T > which fulfills 3N Fadddeev equation 

\T >= > +tPGo\T > , (1) 

where P is defined in terms of transposition operators of three nucleons, P = P12P23 + -P13-P23) Gq 
is the free 3N propagator, and t is the two-nucleon off-shell t-matrix. Knowing \T > the breakup 
as well as the elastic nd scattering amplitudes can be gained by quadratures in the standard 
manner Namely, the transition amplitude for the elastic scattering, < <I>'^ ]^|C/|<l>i^i >, is given 
by [l|,B 

(i>;,iic/|ci>i,i> = (<i>;,i|PGo>i,i> + (ci>;,ilp|r> , (2) 

and for the breakup, < $o|^o|*^*i,i >) by 

($o|f/o 1-^1,1) = {^o\{l + P)\T) . (3) 

The state \^q) = "^(1 ~ P23) \pq) corresponds to a kinematically complete configuration of the 
breakup described by standard Jacobi momenta p and q, and {^'n > is the outgoing state of the 
elastic scattering with changed direction of the relative neutron-deuteron momentum ^0 ' but with 
the same magnitude as in the entrance channel 'I = l^ol- 

Introducing the momentum space 3N partial wave basis \pqc( >= 
\pq{ls)j{Xl/2)I{jI)J{tl/2)T > with the two-body subsystem angular momenta, spin and 
isospin {ls)j and t, coupled together with the corresponding quantum numbers of the spectator 
nucleon (Al/2)/ and 1/2 to the total angular momentum J and isospin T of the 3N system, 
and projecting Eq. ([1]) on these states, we get the system of coupled integral equations in 
two continuous variables p and q. For details of the numerical treatment of that system, and 
particularly of the kernel part {pqa\tPGo\T) , we refer to 
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The 2-nucleon t-matrix conserves the spectator momentum q and all discrete quantum numbers 
except the orbital angular momentum I: 



<pqa\t\p'q'a'> = ^-^^^tl2if-{pp'-E{q) = E-^q^) 



and has a pole in channels a for which two-nucleon subsystem has bound state. 

In the channels \a >= > which contain the 2-body ^Si —^Di states we extract the deuteron 
pole. Thus we define 



.s^jc^t^l /. p. ^^ ^ i\2lf"{p,p'-E{q)) 



for the deuteron quantum numbers Sq, = = 1, = 0, la, la = 0, 2 and keep t as it is otherwise. 
That pole property obviously carries over to the T-amplitude and we define just for the \ a >= {aa > 
channels 

(p^-m = ^,.'"^f> , . (6) 

E + ie- ^q^ - ed 

Since the energy E of the 3N system is determined by the incoming neutron energy Ecm- E = 
Ecm + ed= -^ql + td, the deuteron pole occurs at g = go- 

When beside the deuteron an additional bound state exists in some 2-nucleon partial wave state 
one needs to extract in channels |q > which contain that 2-nucleon state the corresponding pole 
of the t-matrix by performing the same procedure as for the deuteron. Let us assume that this 
state is a bound state of two neutrons in the "^Sq state with the wave function and the binding 
energy e^n, and let us denote by |<I>i,2 >= Wq, <Pnn > the two body channel build on such dineutron, 
from which or to which different reactions can be initiated. 

In the channels \a >= Ich^q > which contain the ^5*0 dineutron we define 

,s^J.t.f / p. - irjf'-iP^P'-^Eiq)) 

tj^ {P^P ■,E{q)) = ——, — = 3 57—— 7 

E + ie-^q^-enn 3^(90-^)+^^ 



and the dineutron pole occurs at g = go = y 9o + ^(^d ~ Cnn)- Again that pole property carries 
over to the T-amplitude and we define for the \a >= Ickisq > channels the amplitude {pqa\T) 
similarily to Eq. ([6]). The numerical treatment of that new pole follows the treatment of the 
deuteron pole [1] and it requires the set of q-points which, in addition to g = go needed for 
numerical treatment of the deuteron pole, contains also q = % point. Since the dineutron occurs in 
the neutron-neutron ^5*0 state it implies charge independence breaking and the resulting difference 
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between ^5*0 nn and np interactions causes that the proper treatment of the nd reactions recjuires 
inclusion of the total 3N system isospin component T = 3/2 for channels a containing ^Sq [9]. 

The existence of ^Sq dineutron increases number of possible reactions with three nucleons what 
in consequence makes that the unitary relation have to be generalized to include those additional 
processes. It has the form 

< '^l,a\U\^l,a' >* - < ^l,a'\U\^l,a > 

= 27ri ^ / d\ < $,-,6|f/|a>i,,, >* 5(4- i^g-) < ^g-b\U\<^l,a > 
6=1,2 

+ 27ri/6 J d'^pd'^q < $o|f/o|$i,a' >* S{Epg - E^) < ^o\Uo\^i,a > (8) 

with a = 1 and 2 for the deuteron and dineutron channels, respectively. One can choose a' = 1 
or a' = 2 and a = 1 or a = 2. For a = a' = 1 this leads on the left side to the forward scattering 
amplitude and on the right to the total cross section. The energies — Eq~^ are given by the 
binding energies of the deuteron Ei = or dineutron E2 = enn- 

The angular distribution for the process n + d — )■ p + dineutron is given by the transition 
amplitude < <l>i,2|f^|<^*i,i > 

^(^n + d^ p + dineutron) = A\2n)^^ V | (<I>i.2| f/ |<i>i,i> |^ (9) 

us i o On 

mpirinm^ 



where PGq and PT contributions to U are given by 

(^'i,2| -PGq^ 1^1,1) = {4>nn,'mp,qo\PGQ^ \(j)d,mn,md,qo\\z) 



2 r 1 .1 2 2 ^ -xn.l 1 ,11 x.l 1 , 1 1 x 



4tt mM'" ' ^"'^^22 ' 2 2 ''22 ' 2 2 

4>nn{\qo + -^qol) ^ {llMmd + mn- nip, -m„ + nip,md) 

1=0,2 

11 , ,11 , 

(--1| - mn,rnp, -rrin + mp)(--0|m„, -m„,0) 

(l>f{\\% + lo|)5^,md+™„-mj,(^gO + lo) (10) 



and 



(«>l,2|P|r) = {<Pnn,mp,qo\P\T) 

= 5Z X] ^hjKo^soQSj^Q{\Q-I\M - ^x' ,^i\M){1-Tq\ - l,-,--)5^Ao,M-M'(9o) 



°° ''^ Gao,a'{qo ,q',x) 

J -I ttS: 



q'^dq' I dx 0n.(vrO (vr2, g', a'|r> , (11) 



'1 "2 



with 



^2 = ]lqi + \q'^ + q'qox • (12) 

It was assumed that the relative neutron-deuteron momentum % in the incoming channel is directed 
along the z-axis. The convention for isospin projections is that for the neutron it is — ^ while for 
the proton In Eq. (jlip channels oq contain the dineutron two-nucleon subsystem quantum 
numbers with isospin to = 1 and its projection = — 1 and the total isospin Tq of the 3N system 
for these channels is To = ^ or To = |. The geometrical coefficient Gao,a'iqo,q' ,x) stems from the 
matrix elements of the permutation operator P 



III. RESULTS 



In the following we will present and compare to the available nd data the theoretical predictions 
for cross sections in elastic nd scattering and breakup assuming different ^Sq nn force. We take the 
CD Bonn [3] potential as the NN interaction and multiplying its ^Sq nn component by a factor A 
generate a number of ^^o nn forces among which some provide binding of two neutrons. In Table |I] 
we show values of the nn scattering length the effective range parameter r^jf and the dineutron 
binding energy a number of A values. Changing A from 0.9 to 1.5 leads to nn *S*o force with 

different, negative as well as positive, values of the scattering length. In order to see if conclusions 
depend from a particular ^Sq nn potential used and from the method applied to generate the nn 
bound state, we performed also calculations with a chiral NN potential in next-to-leading-order 
(NLO) of chiral expansion adjusting its ^Sq nn low energy constants to get a dineutron with 
given binding energy. 



A. Total cross sections 



The results for the nd total cross sections are shown in Fig. [T] and, for a number of energies, 
they are also presented in Table HIl The theoretical predictions obtained with different nn ^Sq 
forces are compared to numerous data taken over many years. Up to ~ 100 MeV there is a 
nice agreement between all data, especially very precise one of Ref. [10], and theory based on the 
CD Bonn potential. When instead of the original CD Bonn ^Sq nn force the modified interaction 
with factor A = 0.9 is taken the resulting cross sections seem to be not excluded by the total cross 
section data. For A = 1.21, with the dineutron binding energy €nn = —144 keV, the predicted total 
cross sections for energies up to w 10 MeV differ from the data by about three standard deviations. 
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At higher energies they clearly lie outside three standard deviations from the data. Increasing the 
factor A to 1.3 or 1.4 leads to total cross section values strongly overestimating the data. 

In Figs. [2] and [3] we compare theoretical predictions for the total elastic scattering and breakup 
cross sections, respectively, with the corresponding data. For the elastic scattering component of 
the total cross section (see Fig. [2]), at energies up to about En ~ 20 MeV, the theoretical predictions 
with different nn ^Sq forces are close to each other and they agree with the data. At energies above 
En ~ 20 MeV cross sections for A > 1 start to deviate from standard CD Bonn and A = 0.9 values 
and the data seems to prefer larger values of A. 

For the total breakup cross sections (see Fig. [3|) the data seems to be compatible with all 
theoretical predicitons with exception of data from Ref. jisl . That data set taken in the region of 
energies 12 MeV < En < 22 MeV clearly advocates the CD Bonn potential predicitons However, 
it does not exclude definitly values of cross sections obtained with A = 1.21. 

At low energies the nd interaction is parametrized by the doublet, ^flnd, and quartet, ^a^^, 
scattering lengths. While '^and is strongly influenced by a 3NF the ^and is practically insensitive to 
such an interaction [3]. In Table Hill we show how these scattering lengths change with modification 
of the ^5*0 nn CD Bonn potential. While the doublet scattering length drastically changes with A 
the quartet scattering length practicaly remains constant under such modifications of the ^Sq nn 
force staying close to the experimental value ^a^^ = (6.35 ± 0.02) fm [17 [. 



B. Elastic neutron-deuteron scattering 

The nd elastic scattering angular distributions are shown in Fig. HI At cm. scattering angles 
©cm. > 45° different theories practically overlapp and agree with the data for all four energies 
shown. Such behaviour is not suprising since at backward angles the exchange term PGq^ , given by 
the deuteron wave function, dominates the elastic scattering transition amplitude. The properties 
of the nn interaction should play decisive role at forward angles. Indeed, at forward angles 
below Qc.m. < 45° differences between theoretical predictions based on different nn ^5o forces start 
to appear and they increase with deacreasing angle. However, in the forward angular region the 
nd elastic scattering cross section data are lacking with exception of En = 14.1 MeV where 5 
data points fall into that region of angles. While two data points at smallest angles support the 
CD Bonn cross sections three other at greater angles prefer the larger values of A. The precise nd 
elastic scattering data at forward angles are required to decide if stronger nn ^5*0 force is allowed. 
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C. Breakup 



Among numerous kinematically complete breakup configurations the largest discrepancies be- 
tween theory and data have been found for the nn QFS and SST geometries. For these configu- 
rations the theoretical cross sections are insensitive to the underlying dynamics and they do not 
change when applying different realistic NN potentials and combining them with available 3NF's. 
Also when instead of the nn QFS one compares the theory with the only one available np QFS data 
set of the nd breakup 2^ a nice agreement is found. QFS and SST configurations are dominated 
at low energies by the ^Sq and ^Si NN force components {g}] which practically saturate the QFS 
and SST cross section at low energies Q]- It would suggests that it is the nn ^Sq force which is 
probably responsible for large discrepancy between data and theory. 

For QFS configurations we show in Fig. [5] the sensitivity of the nn and np QFS configurations to 
the underlying nn ^5*0 force. As expected, changes of that force cause drastic modifications of the 
nn QFS cross section leaving the np QFS practically without modifications. As was shown in 3] 
responsible for such drastic modifications of the nn QFS cross section are changes of the effective 
range parameter induced by factor A. The changes in the nn scattering length practically leave the 
nn QFS cross sections without modifications. 

In contrast to the nn QFS the SST geometry is more stable against changes of the ^5*0 nn 
force. As shown in Fig. [6] changing the factor A does not bring theory closer to the data. While 
A = 0.9 provides smaller SST cross sections than the CD Bonn potential, taking factor A > 1 
and increasing it so that dineutron is formed, leads to cross sections which again are below the 
CD Bonn potential predicitons. Therefore by modifications of the ^Sq nn force it is not possible 
to explain the large discrepancy for SST. Since it is unprobable that the deuteron properties are 
so badly known that the ^Si Di NN force component would require modification, the source 
for that disagreement must be sought elsewhere. One possibility could be the indirect infiuence by 
the dineutron some breakup configurations by contributing in specific regions of a phase-space to 
the breakup background. 

For the ^Sq nn force which allows a dineutron the nn scattering length becomes positive. It 
should have drastic infiuence on the nn final-state-interaction (FSI) of the nd breakup, where the 
two outgoing neutrons, having the same momenta, are strongly interacting in the ^5*0 state. We 
show in Fig.[7]the changes in the FSI peak when the nn scattering length a„„ changes from negative 
to positive values. For the same magnitude of a„„ the nn FSI cross section is strongly diminished 
for the positive sign of ann- The question arises if the existing nn FSI cross section data can be 
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understood when dineutron exists ? 

To answer that question we show in Fig. [5] cross sections for 4 kinematically complete nn FSI 
configurations for which data have been taken and analyzed in Ref. 26|] with the aim to extract the 
neutron-neutron scattering length. The consistent values of a„„ have been found in each of those 
4 configurations with the average value of a^n = 18.7 it 0.7 fm. As can be seen in Fig. [51 again 
changing a„„ to positive values diminishes significantly the nn FSI cross sections. Comparing 
cross sections obtained with A = 1.19 and A = 1.21 to the CD Bonn potential values clearly 
demonstrates that no theoretical analysis of [26^ data, when performed with positive values of 
o-nn, would provide consistent values for the nn scattering length in those 4 geometries. While the 
analysis of 61 = 62 = 43° configuration would probably provide a„n = +21.69 fm, a similar analysis 
of configurations at smaller 9i = 62 would provide distinctly larger positive nn scattering lengths. 

In Fig. [9] we show further 3 FSI configurations for which data have been taken. For d{n, nn)p 



complete breakup the data of Ref. [2j] support the CD Bonn potential cross section predictions. 
Each of 2 complete configurations in that reaction shown in Fig. [9] contain two np FSI peaks. 
The theoretical analysis of these np FSI peaks, if performed with positive values of Onn, would 
provide different values for the neutron-proton scattering length a„p, which in addition would be 
inconsistent with well known a^p experimental value. 

In Fig. Owe show also the configuration for d{n,np)n breakup in which data have been taken 
and analyzed in Ref. [27]. That geometry contains both np and nn FSI peaks. Again, the analysis 
of the np FSI peak, if performed with positive a„„, would provide too large magnitude for a^p. 

To see how our conclusions depend on the NN potential used and on the method applied to 
modify the ^Sq nn force we present in Fig. [8] also cross sections obtained with next-to- leading-order 
(NLO) chiral perturbation theory potential of Ref. [5], including in calculations all np and nn forces 
up to the total angular momentum jmax = 3 in the two-nucleon subsystem. The component 
of that interaction is composed of the one- and two-pion exchange terms and contact interactions 
parametrized by two parameters Ci and Ci 5^ 

veso) = as, + as,{p'+p'') ■ (13) 

Standard values are Cis^ = -0.1557374 x 10^ GeV^ and Cig^ = 1.5075220 x 10^ GeV^ for cutoff 
combinations {A, A} = {450 MeV, 500 MeV} [5]. By multiplying Cig^ by a factor C2CS0) and 
Ci5p by a factor Ci(^S'o), one can induce changes of the nn ^Sq interaction. In Fig. [8]we show two 
predictions based on the NLO potential with negative (a„„ = —17.6 fm - the (magenta) dashed- 
dotted line) and positive (a^^ = +17.5 fm - the (green) double-dashed-dotted line) values of the 
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neutron-neutron scattering length. Comparing them with different CD Bonn potential predictions 
and taking into account differences between their a„„ values it is clearly seen that both potentials 
and methods of changing ^Sq nn interaction lead to the same conclusions. 

FSI region can also be investigated in the uncomplete breakup measurement, in which spectrum 
of the outgoing proton is measured at given lab. angle. In Fig. [TU] we show modifications of the 
outgoing proton spectrum for 14 MeV nd breakup at proton lab. angle 6 = 4:°. Here changing the 
sign of leads to disappearing of the FSI peak. In addition , at lower energies of the outgoing 
proton the modification of ^Sq nn force by factors A > 1 significantly increases the uncomplete 
breakup cross section. 

The analysis of existing nd uncomplete breakup spectra performed in 



29|] indicated on the 



inconsistencies in the experimental uncomplete nd breakup data and revealed unexplained differ- 
ences of more than 25% in regions of the outgoing proton energies where large number of different 
three-nucleon configurations contribute to the cross section. The question arises if existence of the 
dineutron and corresponding modification of the ^Sq nn force can account for that and if the clear 
FSI peaks appearing in the experimental outgoing proton spectra provide evidence for existence 
of dineutron. In order to answer that question a theoretical Monte Carlo analysis of experimental 
spectra, which would provide the angular distribution for the dineutron cross section, is required. 
The resulting angular distribution should then be compared to the theoretical angular distribution 
for n + d ^ p + dineutron transition. However, in view of the results presented above for the com- 
plete nn FSI configurations for which data are available, it seems highly unprobable that analysis 
of incomplete spectra will provide a clear signal for existence of dineutron. 



D. Transition from the neutron-deuteron to the proton-dineutron channel 

For values of factor A = 1.21, 1.3 and 1.4, which allow for the bound ^5*0 state of two neutrons, 
the transition to the proton-dineutron channel is possible. In Fig. [TT]we show angular distributions 
for n+d — )• p+dineutron reaction. The cross sections for that reaction are by an order of magnitude 
smaller than for the nd elastic scattering, with largest cross sections at backward cm. angles for 
low incoming neutron energies. 

In view of the discrepancies found in the nd breakup reaction, especially in SST configuration, it 
is interesting to consider in which phase-space region that hypothetical dineutron state could mostly 
affect the breakup configurations by contributing in an uncontrolled manner to the background. 
To answer that question again Monte Carlo simulations of experimental conditions are required. 
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IV. SUMMARY AND CONCLUSIONS 



We have investigated how far available nd data allow for hypothetical bound state of two 
neutrons and if such dineutron can help to explain the discrepancies between theory and data found 
in some complete nd breakup configurations. To this aim we extended our numerical momentum 
space treatment of 3N Faddeev equations to include in addition to the deuteron also ^Sq dineutron. 
Solution of these equations with modified nn ^Sq CD Bonn force provided predictions for cross 
section in different nd reactions. 

We found that available nd data for the total nd interaction cross section are incompatible with 
the existence of a dineutron with binding energy of absolute value greater than 100 keV. The 
data for the total elastic scattering and breakup cross sections do not exclude such a possibilty. 
Also data for the nd elastic scattering angular distribution cannot decisively exclude such a state. 
However, in this case the precise data taken at forward angles, if available, could provide more 
constraints on existence of a dineutron. 

The modifications of the ^5*0 nn force component cannot provide explanation for the drastic 
discrepancy between theory and data for the SST geometry of the nd breakup. Allowing for 
dineutron provides even smaller SST cross sections, increasing thus that discrepancy. 

The transition from negative to positive nn scattering lengths leads to drastic modifications of 
the FSI cross sections. In the outgoing proton spectra of the uncomplete nd breakup the positive 
scattering length leads to strong diminishing of the FSI peak at maximal proton energies. The 
carefuU Monte Carlo theoretical analysis of existing proton spectra is required to get answer if these 
spectra provide a clear signal for the existence of dineutron. However, complete FSI configurations 
for which data have been taken exclude positive values for a„„. 
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TABLE I: The dineutron binding energy e„„, the nn scattering length a„„ and the effective range parameter 
Veff for different factors A by which the nn ^Sq component of the CD Bonn potential was multiplied. 



\ 

A 


enn [MeVJ 


Unn [imj 


^ f-P, 

^eff [11 


0.9 




-8.25 


3.12 


1.0 




-18.80 


2.82 


1.19 


-0.099 


-1-21.69 


2.39 


1.21 


-0.144 


-1-18.22 


2.35 


1.3 


-0.441 


-1-10.95 


2.20 


1.4 


-0.939 


+7.87 


2.07 
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TABLE II: The theoretical (evaluated at the nucleon laboratory energy Eth) and experimental (taken at 
Eexp) nd total cross sections. Theoretical values were obtained with the CD Bonn potential the nn 
component of which was multiplied by a factor A. 



Eth 
[MeV] 


^exp 

[mb] 


[MeV] 


[mb] 


A=1.0 
"th 

[mb] 


A=1.21 

[mb] 


A=1.3 
"th 

[mb] 


A=1.4 

[mb] 


8.0 


1207 ± 13 
1213.3 ± 5.58 
1224 ± 10 


8.0 [32] 
8.038 [30] 
8.0 [nj 


1203.4 


1205.6 


1258.5 


1301.4 


1353.5 


10.0 


1055 ± 10 
1051.1 ± 6.9 
1045.0 ± 3.4127 


10.0 [32] 
10.026 [30] 
9.9218 [10] 


1026.4 


1036.1 


1089.9 


1123.7 


1162.5 


13.0 


867 ± 12 


12.995 [32] 


837.96 


851.76 


900.91 


926.56 


954.72 


14.1 


803 ± 14 
790 ± 20 

809 ± 6 
778 ± 22 
806 ± 6 

810 ± 30 


14.1 [34] 
14.1 [33] 
14.1 [35] 
14.1 [36] 

14.1 [36] 

14.2 [37] 


783.94 


798.25 


845.37 


868.91 


894.52 


19.0 


627.96 ± 12.16 
632 ± 14 


18.932 [30] 
19.01 [31] 


603.47 


617.55 


655.92 


673.20 


691.76 


26.0 


455 ± 12 
451.47 ± 17.72 


26.015 [32] 
26.082 [30] 


444.41 


456.18 


485.43 


497.83 


511.16 


42.5 


267.7 ± 3.9 


42.5 [38] 


259.24 


266.35 


283.32 


290.27 


297.88 


65.0 


166.5 ± 2.9 
161.7 ± 2.8 
168.27.0 ± 0.48333 


63.5 [38] 
66.5 [38] 
65.039 [1^ 


157.24 


160.95 


170.27 


173.96 


178.13 
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TABLE III: The doublet '^Und and quartet *and nd scattering lengths for different factors A by which the nn 
^5*0 component of the CD Bonn potential was multiplied. The calculations have been done with all partial 
waves with 2N total angular momenta up to jmax = 3 included. 



A 



0.9 

1.0 

1.21 

1.3 

1.4 



1.51485 
0.93174 
-0.43567 
-1.18887 
-2.37605 



6.34602 
6.34600 
6.34596 
6.34593 
6.34589 
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FIG. 1: (color online) The total cross section for the neutron-deuteron interaction as a function of the 
neutron lab. energy. Different lines show sensitivity of the total cross section to the changes of the nn 
^5*0 force component. Those changes were induced by multiplying the ^5*0 nn matrix element of the CD 
Bonn potential by a factor A. The solid (blue) line is the full result based on the original CD Bonn potential 
(A = 1.0) and all partial waves with 2N total angular momenta up to jmax = 3 included. The (black) dotted, 
(red) short-dashed, (maroon) long-dashed, and (green) dashed-double-dotted lines correspond to A = 0.9, 
1.21, 1.3, and 1.4, respectively. The (magenta) circles, (red) x-es, and (maroon) diamonds are nd data of 
Ref. Q, [ul, and 12 1, respectively. 
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FIG. 2: (color online) The total neutron-deuteron elastic scattering cross section as a function of the neutron 
lab. energy. Different lines show sensitivity of that cross section to the changes of the nn ^5*0 force component. 
For their description see FigHJ The (maroon) diamonds are nd data of Ref. . 
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FIG. 3: (color online) The total neutron-deuteron breakup cross section as a function of the neutron lab. 
energy. Different lines show sensitivity of that cross section to the changes of the nn ^5*0 force component. 
For their description see Fig |Tl The (red) crosses, (green) circles, (blue) squares, and (maroon) diamonds 
Q, ll5| and respectively. 
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FIG. 4: (color online) The neutron-deuteron elastic scattering angular distributions da/dQ at a number of 
incoming neutron lab. energies. Different lines show sensitivity to the changes of the nn ^5*0 force component. 
For their description see FigH] At _E„ = 8, 10 and 14 MeV the (magenta) circles and (red) x-es are nd data 
of Ref. Q and 0, respectively. At En = 14 MeV the (maroon) stars, (green) triangle-down, and (black) 
triangle-up are nd data of Ref. H, Q, and 0, 

respectively. At £"„ = 65 MeV the (blue) circles are 



En = 66 MeV nd data of Ref. 
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FIG. 5: (color online) The cross section (f' a / (KlidVl2dS as a function of the S-curve arc-length in the 
j^iab _ 26 MeV nd breakup reaction d[n^ nn)p for the QFS nn (left) and np (right) kinematically complete 
configurations of Ref. [2^. For description of lines see Fig.[T] 
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FIG. 6: (color online) The cross section a / <KlidVl2dS as a function of the S-curve arc-length in the 
E'^"'' = 13 MeV nd breakup reaction d{n,nn)p for SST configuration with the lab. angles of two detected 
neutrons 9i = 62 ~ 52.8° and 4)12 = 180°. For description of lines see Fig.[T] The (magenta) solid dots and 



(red) x-ses are nd data of Ref. [2, 



and 



251 ■ respectively. 
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FIG. 7: (color online) The cross section a / dVli(Kl2dS for the E]!^^ = 13 MeV nd breakup reaction d{n, np)n 
as a function of the S-curve length for nn FSI configuration. For description of lines see Fig. [T] 
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FIG. 8: (color online) The cross section a / (KlidVl2dS as a function of the S-curve arc-length in the 
— 13 MeV nd breakup reaction c?(n, nn)p for 4 FSI nn geometries. Different lines show sensitivity of 
the cross section to the changes of the nn ^5*0 force component. Those changes were induced for the dotted 
(black), (red) dashed and (blue) solid lines by multiplying the ^Sq nn matrix element of the CD Bonn 
potential by a factor A. The dotted (black) line is the full result based on the original CD Bonn potential 
(A = 1.0) and all partial waves with 2N total angular momenta up to jmax — 3 included. The (red) dashed 
and (blue) solid lines correspond to A = 1.19 and 1.21, respectively. The (magenta) dashed-dotted and 
(green) double-dashed-dotted lines show results of Faddeev calculations based on NLO chiral perturbation 
theory potential and all partial waves with 2N total angular momenta up to jmax = 3 included. They differ in 
the nn ^Sq force which for the (magenta) dashed-dotted line was obtained with the constants Ci(^5o) ~ 1.0 
and C2(^5o) = 1.0 (original NLO potential, see text for explanation) leading to a„„ = —17.6 fm and r^ff = 
2.75 fm. For the (green) double-dashed-dotted hue the constants Ci(^S'o) = 1.50 and C2(^S'o) = 1.29415, 
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FIG. 9: (color online) The cross section a / (KlidVl2dS as a function of the S-curve arc-length in the 
i?^"'' = 13 MeV nd breakup reaction d{n, nn)p for 3 FSI geometries. Different lines show sensitivity of that 
cross section to the changes of the nn ^Sq force component. For their description see FiglHl 
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FIG. 10: (color online) The cross section (Pa/dflidEp for the — 14 MeV uncomplete nd breakup 

reaction d(n,p)nn as a function of the outgoing proton lab. energy at the proton lab. angle Oi = 4°. For 
description of lines see Fig. [TJ 
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FIG. 11: (color online) The angular distributions da/dQ for d{n,p)dineutron reaction at a number of 
incoming neutron lab. energies. The (red) short-dashed, (maroon) long-dashed, and (green) dashed-double- 
dotted lines correspond to the factor A = 1.21, 1.3, and 1.4, respectively. 
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